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We report the results of a search for T and CP violation in B°-B° mixing using an inclusive 
dilepton sample collected by the BABAR experiment at the PEP-II B Factory. The asymmetry 
between £ + £ + and £~£~ events allows us to compare the probabilities for B° — > B° and B° — > B° 
oscillations and thus probe T and CP invariance. Using a sample of 23 million BB pairs, we measure 
a same-sign dilepton asymmetry of A T / CP = (0.5 ± 1.2(stat) ± 1.4(syst)) %. For the modulus of the 
ratio of complex mixing parameters p and q, we obtain \q/p\ = 0.998 ± 0.006(stat) ± 0.007(syst). 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Since the first observation of CP violation in 1964 [Q , investigation of CP violation through the measurements 
the kaon system has provided many other results prob- of the unitarity triangle angles a, (3 and 7, the BABAR 
ing the CPT and T discrete symmetries . Beyond the experiment can investigate T and CP violation purely in 
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mixing. 

The physical states (solutions of the complex effective 
Hamiltonian for the B°-B° system) can be written as 

\Bl H )=p\B )±q\B°) 

where p and q are complex mixing parameters with the 
normalization \p\ 2 + \q\ 2 = 1. 

The CPT invariant asymmetry, A T / CP , between the 
two oscillation probabilities P(B° -> B°) and P{B Q -> 
B°) probes both T and CP symmetries and can be ex- 
pressed in terms of p and q: 

_ P{B° -> B°) - P{B° -> B°) 
T/CP ~ P(B° -> 5°) + P(£° -> £°) 

_ i - k/p| 4 
i + k/p| 4 ' 

Standard Model calculations |3| predict the size of this 
asymmetry to be at or below 10~ 3 . Therefore, a large 
measured value could be an indication of new physics. 

Inclusive dilepton events representing 4% of all 
T(45) — > BB decays provide a very large sample with 
which to study T and CP violation in mixing. The fla- 
vor of each B meson is tagged by the charge of the lep- 
ton. Assuming AB = AQ and CP invariance in the 
direct (b — > £) semileptonic decay process, the asym- 
metry between same-sign lepton pairs, £ + £ + and £~£~, 
allows a comparison of the two oscillation probabilities 
P(B a -v B°) and P(B° -> B°). The asymmetry A T/CP 
for direct same-sign dileptons is time independent. How- 
ever, in this analysis, the time difference At between the 
two B meson decays is used to discriminate the direct lep- 
tons from the cascade leptons produced in (b — > c — > £) 
transitions. 

The measurement of A T iq p reported here is performed 
using events collected by the BABAR detector jj] from 
e + e _ collisions at the PEP-II asymmetric-energy B Fac- 
tory between October 1999 and October 2000. The in- 
tegrated luminosity of this sample is 20.7 fb _1 recorded 
at the T(45) resonance ( "on-resonance" ) and 2.6 fb _1 
recorded about 40 MeV below the T(45) resonance ("off- 
resonance"). BB pairs from the T(4S) decay move 
along the high-energy beam direction (z) with a nomi- 
nal Lorentz boost (fij) — 0.55. 

Lepton candidates must have at least 12 hits in the 
drift chamber (DCH), at least one z-coordinate hit in 
the silicon vertex tracker (SVT), and a momentum in 
the T(4S) center-of-mass system (CMS) between 0.7 and 
2.3 GeV/c. Electrons are selected by requirements on 
the ratio of the energy deposited in the electromagnetic 
calorimeter (EMC) and the momentum measured in the 
DCH, on the lateral shape of the energy deposition in 
the calorimeter, and on the specific ionization density 
measured in the DCH. Muons are identified through the 
energy released in the calorimeter, as well as the strip 



multiplicity, track continuity and penetration depth in 
the instrumented flux return (IFR). Lepton candidates 
are rejected if they are consistent with a kaon or pro- 
ton hypothesis according to the Cherenkov angle mea- 
sured in the detector of internally reflected Cherenkov 
light (DIRC) or to the ionization density measured in 
the DCH. The electron and muon selection efficiencies are 
about 92% and 75%, with pion misidentification proba- 
bilities around 0.2% and 3%, respectively. 

Non-BB events are suppressed by requiring the ra- 
tio of second to zeroth order Fox- Wolfram moments 0] 
to be less than 0.4. In addition, the residual contam- 
ination from radiative Bhabha and two-photon events 
is reduced by requiring the squared invariant mass of 
the event to be greater than 20GeV 2 /c 4 , the event apla- 
narity to be greater than 0.01, and the number of charged 
tracks to be greater than four. Electrons from photon 
conversions are identified and rejected with a negligi- 
ble loss of efficiency for signal events. Leptons from 
J/tp and ip(2S) decays are identified by pairing them 
with other oppositely-charged candidates of the same- 
lcpton species, selected with looser criteria. We re- 
ject the whole event if any combination has an invari- 
ant mass within 3.037 < M (£+£') < 3.137 GeV/c 2 or 
3.646 < M (£+£') < 3.726 GeV/c 2 . 

To minimize the wrong flavor tags due to leptons from 
cascade charm decays, we use a neural network algo- 
rithm that combines five discriminating variables. These 
are calculated in the CMS (see Fig. |l|) and are the mo- 
menta of the two leptons with highest momentum, and 
P2, the total visible energy E to t, the missing momentum 
Pmiss of the event, and the opening angle between the 
leptons, 6i2- The first two variables, pi and p%, are very 
powerful in discriminating between direct and cascade 
leptons. The last variable, #12, efficiently removes direct- 
cascade lepton pairs coming from the same B and further 
rejects photon conversions. Some additional discriminat- 
ing power is also provided by the other two variables. In 
order to be insensitive to the Monte Carlo, the fraction of 
cascade leptons is determined from a fit to the same-sign 
and opposite-sign dilepton data. 

In the inclusive approach used here, the z coordinate of 
the B decay point is the z position of the point of closest 
approach between the lepton candidate and an estimate 
of the T(iS) decay point in the transverse plane. The 
T(45) decay point is obtained by fitting the two lep- 
ton tracks to a common vertex in the transverse plane, 
which is constrained to be consistent with the beam- 
spot position. The time difference, At, between the two 
B meson decays is determined from the absolute value, 
Az, of the difference in z between the two B decays 
by At = Az/(/3j)c. The background events (cascade 
leptons from unmixed B°B° events and B + B~ events, 
and non-Si? events) are most prominent at low Az (see 
Fig. |^). Therefore, a requirement of Az > 200 /im al- 
lows us to eliminate about 50% of background without 
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FIG. 1: Distributions of the discriminating variables (a) p\, 
(b) pi, (c) Etot, (d) Pmiss and (e) 612, for data (dots) and 
Monte Carlo events (histograms). The contributions from 
direct-direct pairs, direct-cascade or cascade-cascade pairs, 
and pairs with one or more fake leptons are shown for the 
Monte Carlo samples. 



dramatically decreasing the signal efficiency. Finally, in 
the measurement of A T j CP , the dilution factor due to 
remaining background is corrected as a function of At. 

Application of the selection criteria described above 
results in a sample of 20,381 same-sign dilepton events, 
consisting of 5,252 electron pairs, 5,152 muon pairs and 
9,977 electron-muon pairs. For Az > 200 /im, the frac- 
tion of now-BB events, measured with the off-resonance 
data, is 4.3% with a charge asymmetry of (—5 ± 10)%; 
the main BB backgrounds, determined from Monte Carlo 
simulation, include 24% of one direct lepton paired with 
a cascade lepton from the other B, 10% of fake leptons 
from the other B, 2% of fake leptons from the same B 
and 2% of leptons from J/0 or resonance decays. 

Since the asymmetry A T j CP is expected to be small, 
we have carefully determined the possible charge asym- 
metries induced by the detection and reconstruction of 
electrons and muons. The three sources of charge asym- 
metry in the selection of lepton candidates come from 
differences, for positive and negative particles, in track- 
ing efficiency £% ack , in particle identification efficiency 
£p id , and in misidentification probability rf^ id . Indepen- 
dent samples are used to estimate these efficiencies and 
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FIG. 2: Distribution of the same-sign dileptons as a function 
of Az; the curve superimposed on the dots is determined from 
a fit to the same-sign and opposite-sign dileptons; the solid 
and dotted lines represent respectively the signal component 
(B°B° or B°B° pairs) and the background component (cas- 
cade leptons from unmixed B°B° and B + B~ events, leptons 
from J/ip , resonance decays and non- BB events). 

probabilities as a function of several charged track pa- 
rameters Xi'. total or transverse momentum, and polar 
and azimuthal angles in the laboratory frame. The num- 
bers of "detected" positive and negative leptons N^ et are 
related to the numbers of true leptons N true by the equa- 
tion 

N deM>P*) =N? rue (x i ,p*) ■ ef rack (xi) ■ [e% d { Xl ) + 

r(ir,P*)-vf id (K,Xi)+r(K,p*)- (2) 
V& d ( K > x i) + r (P,P*) ■ V%(P, Xi)] , 

where r(n,p*), r(K,p*) and r(p,p*) are the relative 
abundances of hadrons (tt, K, and p) with respect to 
the lepton abundance for a given p* (the momentum 
of the track in the CMS). These quantities are ob- 
tained from BB Monte Carlo events, after applying the 
event selection criteria with perfect particle identifica- 
tion. To correct for charge asymmetries in the lepton 
detection, we apply a weight proportional to the ratio 
N true (xi,p*)/N d t et (x i ,p*), for each lepton in the sample. 

Using tracks selected from multi-hadron events, the 
tracking efficiencies £^.^(2^) for positive and negative 
particles are determined by computing the ratio of the 
number of SVT tracks with at least 12 DCH hits as re- 
quired in the dilepton selection, divided by the initial 
number of SVT tracks. These tracking efficiencies are 
tabulated as a function of transverse momentum and po- 
lar and azimuthal angles. The charge asymmetry correc- 
tion is less than 0.1% on average in the relevant momen- 
tum range. 

The identification efficiencies £^ id {xi) are measured as 
a function of total momentum and polar and azimuthal 
angles, with two control samples consisting of ee — > 
eeee (with 77 — > ee) and radiative Bhabha events for elec- 
trons, and with a ee — > eefifj, (with 77 — > fifi) control 
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sample for muons. The misidentification probabilities 
(hadron, Xj) are determined using control samples of 
kaons produced in D* + — > ir + D° — > tt + K~it + decays 
(and charge conjugate), pions produced in ifg — > 7r + 7r~ 
decays, and one-prong and three-prong r decays, and 
protons produced in 7l decays. 

For the electrons, the charge asymmetry in the par- 
ticle identification efficiency reaches (0.5-1.0)% in some 
regions of the lepton phase space. The impact of the 
charge asymmetry in misidentification is negligible be- 
cause the absolute misidentification probability for pions 
is extremely small (~ 0.2%). However, the A control sam- 
ple indicates a very large misidentification probability for 
antiprotons with momentum ~ 1 GeV/ c. Such an effect 
is due to the annihilation of antiprotons with nucleons 
in the calorimeter, which produces a signature similar to 
that of an electron. The impact of this effect is balanced 
by the low relative abundance of antiprotons in B decays. 
Overall, antiprotons induce a charge asymmetry of order 
0.1% and a correction is applied for this effect. 

For the muons, the ee/i/j, control sample shows that 
the charge asymmetry in the efficiency reaches 0.5%. 
The misidentification probability for pions is much larger 
(~ 3%) than in the case of electrons but there is no 
indication of any charge asymmetry induced. On the 
other hand, the kaon misidentification distribution shows 
a charge asymmetry at the level of (10-20)% due to the 
difference between cross sections for K + and K~ meson 
interactions with matter for momenta around 1 GeV/c. 

Equation [l] is applicable for pure signal (direct leptons 
from B°B° and B Q B° events). However, the dilepton 
sample is contaminated by cascade leptons from B + B~ 
and unmixed B°B° events, non-Si? events, and J/0 de- 
cays (see Fig. ||). Assuming no charge asymmetry in 
the background and assuming CP invariance holds in di- 
rect semileptonic B decays, we can write the measured 
asymmetry A™^ s p (see Fig. ||) in terms of the number of 
events N as 

A meas ( ... _ N (£+£+, At) - N £~ , At) 

S(At) (3) 
T/CP S(Ai) + B(Ai) ' 

where S(At) and B(At) are the numbers of sig- 
nal and background events respectively. Thcrc- 

re- 



forc, 



extraction of a value for A 



T/CP 



quires a determination of the dilution factor 
S(At) I \S{At) + B(At)] . The asymmetry between 
same-sign dileptons is corrected for the background 
dilution using the time dependent probability density 
functions shown in Fig. ^. These probability density 
functions are obtained with a fit to data for the same- 
sign and opposite-sign dilepton samples with the value 
of Arrid fixed to the world average value §. This fit 
is similar to that used in the measurement of Amj 



with dilepton events 0: it determines the corrections 
to the resolution function extracted from Monte Carlo 
simulations, the fraction of cascade leptons, the average 
lifetime of the charm component for cascade leptons, the 
fraction of cascade leptons, and the fraction of charged 
B events. In addition, the fraction of non-SB events is 
measured from off-resonance data. From a \ 2 fit to the 
distribution of the asymmetry as a function of At for 
the same-sign dileptons with Az > 200 /im (see Fig. ||), 
we extract A T/CP = (0.5 ± 1.2) %. 
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FIG. 3: Corrected same-sign dilepton asymmetry as a func- 
tion of At. The line shows the result of the fit for the dileptons 
with Az > 200 /im. 

The systematic uncertainties related to the detection 
charge asymmetry both for tracking and lepton identifi- 
cation are determined using direct leptons from semilep- 
tonic B decays. This sample has the same topology and 
kinematics as the leptons from dilepton events. This 
single-lepton charge asymmetry, sensitive to the charge 
asymmetry due to detection bias, may also be affected 
by the real physical asymmetry A T / CP in the dilepton 
events. But, in practice, the effect introduced by A^/cp 
is suppressed by more than one order of magnitude and is 
therefore neglected. With the 1999-2000 data set, we se- 
lect roughly 1.5 million electrons and 1.5 million muons. 
After subtraction of scaled off-resonance data and after 
applying a correction weight derived from Eq. |^, we mea- 
sure the charge asymmetries to be (—0.30 ± 0.14)% for 
the electrons and (—0.35 ±0.17)% for the muons. We as- 
sign these residual asymmetries ±0.30% and ±0.35% as 
systematic errors due to charge asymmetry in detection 
efficiencies. With the dilution factor correction, the total 
systematic errors related to the charge asymmetry in the 
detection are ±0.5% and ±0.6% for electrons and muons, 
respectively. 

The assumption of no charge asymmetry in the back- 
ground is confirmed by the off-resonance data where the 
charge asymmetry ( — 5 ± 10)% is consistent with zero and 
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leads to a ±0.7% uncertainty on the A T / C p measure- 
ment. In addition, the charge asymmetry of the events 
with Az < 100 fim, which contain 85% background (cas- 
cade leptons from i3 ± and unmixed B°), is (1.2 ± 1.4)%, 
also consistent with zero. From this asymmetry, we can 
constrain to ±0.9% the uncertainty on A T / CP due to a 
possible charge asymmetry in the decays producing the 
cascade leptons. If we assume CP invariance in the de- 
cays producing the cascade, this uncertainty vanishes. 

The background dilution correction is measured with 
the data from the full dilepton sample with the value of 
Arrid fixed to the world average value The uncer- 
tainty on the ratio B/S leads to a ±3% multiplicative 
error on Ap/cp, which is negligible. A possible dilution 
of A T /cp due to double mistag is neglected because the 
probability of double mistag is at the level of only 1%. 



10 -3 , which is the most stringent test of T and CP vi- 
olation in B°-B° mixing to date and is consistent with 
previous measurements ||. 
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TABLE I: Summary of systematic uncertainties on A t /cp- 



Type of systematic error 


<t(A t/cp )(%) 


[1 


Electron charge asymmetry in the detection 


0.5 


[2 


Muon charge asymmetry in the detection 


0.6 


Non-SB background charge asymmetry 


0.7 


[3 


BB background charge asymmetry 


0.9 


Correction of the background dilution 


0.01 




Total 


1.4 





In conclusion, we measure A T / C p = (0.5 ± 1.2(stat) ± 
1.4(syst)) % where the total systematic uncertainty is the 
quadratic sum of the systematic uncertainties listed in 
Table |. From Eq. [I], the A T j CP asymmetry gives the 
modulus of the ratio of complex mixing parameters p 
and q equal to 

\q/p\ = 0.998 ± 0.006(stat) ± 0.007(syst). 

This measurement can be translated into a measurement 
of the CP violating parameter ep = (p — q)/{p + q)- We 
obtain Re(e B )/(l + |e B | 2 ) = (1.2±2.9(stat)±3.6(syst)) x 
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